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Abstract 

►^ I These notes are the second part of a common course on Renormalization Theory 
fK^ ■ given with Professor P. da VeigafJ. I emphasize here the rigorous non-perturbative or 
fSJ I constructive aspects of the theory. The usual formalism for the renormalization group in 
^D i field theory or statistical mechanics is reviewed, together with its limits. The constructive 
CN I formalism is introduced step by step. Taylor forest formulas allow to perform easily the 
^r^ I cluster and Mayer expansions which are needed for a single step of the renormalization 
Qv I group in the case of Bosonic theories. The iteration of this single step leads to further 
■^ . difficulties whose solution is briefly sketched. The second part of the course is devoted to 
'^ I Fermionic models. These models are easier to treat on the constructive level, so they are 
I . very well suited to beginners in constructive theory. It is shown how the Taylor forest 
i-G ' formulas allow to reorganize perturbation theory nicely in order to construct the Gross- 
ed . Neveu2 model without any need for cluster or Mayer expansions. Finally applications of 
H I this technique to condensed matter and renormalization group around Fermi surface are 
• • . briefly reviewed. 

'V^ ■ 1 The Renormalization Group: an overview 



c^ 



1.1 Functional Integration and its problems 

In this section we restrict ourselves to the bosonic (p'^ field theory in d Euclidean space time 
dimensions. The model, introduced in P. Da Veiga's lectures, is defined by the (formal) 
measure 

d^c('A)e"^^*\ S{ct)) = \ I (t)'^{x)(fx (1.1) 



where d^c{4') represents the Gaussian measure for the free field. Gaussian measures are 
characterized by their covariance, or propagator, which for a massive theory is, in Fourier 
space: 

C{p) = {p^ + mY\ (1.2) 

and S is the (bare) action. In dimension d = 2, 3 the model is superrenormalizable, and 
its rigorous construction was the first major achievement of constructive theory [GJ]. In 
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dimension d = 4 the model is just renormalizable, hence a more general action, including 
mass and wave function counterterms is needed: 



5(0) = A / 0*(x)d'*x + fi (l)^{x)(fx + a {V(l)Y{x)(fx (1.3) 

We recall the three main problems associated with giving a mathematical precise sense to 
the formal measure (1.1). 

A) The ultra-violet problem The propagator (1.2) has not a well defined kernel in direct 
space for d > 2, since J d'^p{p'^ + m?)~^ is not absolutely convergent. This does not prevent 
mathematically to define the Gaussian measure associated to this propagator, for instance 
through Minlos theorem [GJ], since the propagator is positive in Fourier space. However it 
entails that this measure is supported by distributions. The immediate consequence is that 
a local interaction like (p^{x) is ill defined, since multiplication of distributions in general 
is ill-defined. Therefore we must truncate the propagator at high momenta by means of a 
cutoff, introducing for instance 

C,(p)= re-"(p'+"^')da (1.4) 

J K 

which has a well defined kernel 

and to let K ^ later. 

As a consequence the support of duc^i^P) is now made of perfectly smooth functions, and 
0^(x) becomes well-defined. An other popular ultraviolet cutoff, especially for gauge theories, 
is the lattice cutoff, but I will not use it here. 

B) The infra-red problem Even with an ultraviolet cutoff, if 0^(x) is well defined, with 
probability one it is not decreasing at infinity, hence it is certainly not integrable on IR . 
Hence the quantity S{(j)) = /r^ (f)'^{x)d'^x with probability 1 is ill-defined. The solution is to 
restrict the interaction to a finite volume (usually a box A of size L), and to let L ^ cx) later. 

The measure with these two cutoffs becomes 

dficMe'^^^^'^"^'"'" (1.6) 

C) The large field problem In (1.6) the integrand is now well defined on a set of measure 
one. But it does not mean that it can always be integrated! Even in one-dimensional 
integration, F{X) = J_^ e~^ /2-Xx ^^ converges only for A > 0. We cannot hope the infinite 
dimensional integral (1.6) to be better behaved than this one dimensional integral. Therefore 
Bosonic functional integrals require some stability for the potential (here e.g. A > 0). As 
dicussed by P. da Veiga, the perturbation series for F do not converge. F is not analytic, 
but Borel summable. This is also the best we can hope for the 0^ theory, and what has been 
proved in dimension 2 or 3. 

Convergence of the functional integral itself, and the divergence of perturbation theory 
can be considered as generic "large field" problems, because they are related to the fact that a 



bosonic field is an unbounded variable. Physically a large field corresponds to a large number 
of excitations or particles being produced, and large field problems are generic in Bosonic 
theories because Bosons, instead of Fermions, can pile up in large numbers at the same place. 

1.2 Thermodynamic limit 

Having reviewed the main mathematical problems of field theory in the functional framework, 
we remark that cutoffs by themselves alone cannot solve any problem. For instance the 
theory with cutoffs does not satisfy any reasonable axioms. Some manipulations have to be 
performed, cancellations have to occur, so that the limits k ^ and L ^ oo which looked 
hopeless at first sight become finite and well defined. 

The conceptually simplest of these manipulations is the thermodynamic limit A ^ IR , 
which is particularly easy when the theory has a mass. We know from statistical mechanics 
that only the "intensive" or thermodynamic quantities should have a limit as A ^ cx), the 
extensive ones, proportional to the volume, should diverge. Hence the only "manipulation" 
in this case is to restrict our attention to quantities such as the pressure p or the normalized 
Schwinger functions Sn of the theory: 

p= hm ^logZ(A); Z(A) = /e-^^Wd^c(0) (1-7) 

A^oo |A| J 

Sn = lim -i- /,/.(xi)...(/)(x^)e-^^WdMc(0) (1-8) 

In these quotients both quantities should diverge as A ^ co, but their ratio should have a 
finite limit. This is clear in perturbation theory: the series for extensive or unnormalized 
quantities correspond to general graphs, which may contain in particular vacuum graphs 
which by translation invariance, give rise to infinite integrals when the spatial integration is 
over all IR . The intensive quantities such as p or Sn have power series made of connected 
graphs, with at least one vertex fixed, so that if the propagator itself is integrable, as is 
the case in a massive theory, the corresponding amplitudes are finite in the thermodynamic 
limit. In the next sections we will see how the thermodynamic limit can be performed 
non-perturbatively. For Bosons it requires a cluster and a Mayer expansion, whereas for 
Fermions, it requires almost nothing: simply reorganizing the finite dimensional integrals of 
perturbation series in terms of trees rather than Feynman graphs. 

1.3 Renormalization 

The ultraviolet problem increases in difficulty when d increases. In general we cannot expect 
a finite k ^ limit for quantities such as the Schwinger functions if we keep the bare 
parameters fixed as k ^ 0. But physically these bare parameters are not observable in low 
energy realistic accelerators! P. da Veiga has already explained, at least at the perturbative 
level, how to get rid of the ultraviolet divergences that occur in the Feynman graphs of the 
theory. For d = 4 where the theory is just renormalizable, perturbation theory suggests that 
if we allow the bare parameters of the theory A, fi and a in (1.3) to depend on the cutoff 
K in such a way as to diverge as k ^ 0, cancellations can take place so that the Schwinger 



functions Sn might have a well defined limit as k ^ 0. I say that it only suggests this, since 
bare quantities in standard perturbative renormalization are only defined as formal power 
series in the renormalized coupling. These series are not necessarily summable, as explained 
in da Veiga's lectures, both because of the large number of graphs, and because perturbative 
renormalization creates additional sources for divergences, namely renormalons. It is very 
risky to speculate whether or not an infinite non summable series of diverging terms really 
produces a quantity diverging with the cutoffs (in fact in asymptotically free cases it does 
not!). 

I will not add further remarks to the lectures of da Veiga on perturbative renormaliza- 
tion, except to recall that historically many mathematicians felt renormalized perturbation 
theory "pulls the infinities of quantum field theory under the rug" by hiding them into the 
bare constants. Physicists also felt in the early days of renormalization theory that it was not 
"natural" ; in particular it did not explain why nature seems to prefer renormalizable field the- 
ories. Professor Dyson, one of the founding fathers of renormalization, once told me that for 
him the main surprise of the theory was that it lived so long! In short perturbative renormal- 
ization, although beautiful and encouraging, and certainly useful for many phenomenological 
computations (QED, QCD in some regimes...) is a partial and unsatisfactory answer to the 
problems of quantum field theory both for mathematicians and for physicists. 

1.4 The Renormalization Group 

A key progress in the theory of renormalization was the introduction of the philosophy of 
the Renormalization Group (hereafter called RG) by Wilson and followers. The basic idea 
is the following: since the limit «; ^ is so hard to grasp, let us not perform it at once, 
but in steps. A single step can be well understood, and the problem is reduced to the 
hopefully simpler problem of iterating many times a relatively simple transformation. Hence 
RG does not solve renormalization, but it replaces it by a problem in dynamical systems. This 
seems an almost trivial idea, but it is really also a tremendous change in perspective which 
immediately led to progress. For instance it replaced the old "Landau ghost" or its modern 
version, the renormalon problem, by a different question: is the theory asymptotically free 
(or safe), so that the flow for the coupling constant remains bounded at all scales? Also the 
renormalization group philosophy can be applied directly to infrared problems and statistical 
mechanics with a lot of success. 

As you surely know, the ^l theory for the right sign of the coupling constant is not 
asymptotically free. At least until now this means that we do not know how to perform 
its limit K ^ and end up with a non trivial interacting theory (although (/)f remains 
interesting for pedagogical discussions). But the good news are that other theories such 
as gauge theories or the two dimensional Gross-Neveu and a models are asymptotically free. 
Also the renormalization group explained quickly why nature seemed to prefer renormalizable 
theories: indeed in a generic interaction at a very high momentum scale, non renormalizable 
terms are irrelevant; washed out by the RG flow, their presence cannot be detected in the 
effective theory at low energy. 

Let us briefly sketch what would be the RG strategy for the (pf theory: 

The key idea is to split the space of all frequencies into discrete slices, following a geometric 



progression. This can be done on lattices with the popular tool called block-spinning, but we 
can also simply cut the propagator into "momentum slices". Taking some fixed ratio M for 
the geometric progression, and k = M^'^f, p being an integer, we have: 

cM = t.c'{p) (1-9) 

j=0 

C^{p)= e-°(p'+™')da ; C^{p)= e-"(p'+"^')(ia for j > 1 (1.10) 

There is a corresponding separation of the field into a sum of independent random variables: 
(j) = '^j<pj, <pj distributed according to C^ . (j)j is called the (j'-th) fluctuation field and 
J2k=o 'Pk the (j'-th) background field. 
Exercise 1.1 Prove that for some constant K: 

C^{x,y) < KM^^e-^^''^-^'/-^ (1.11) 

Now introduce the operation * by 

fi,*Z{cl>) = Jdfic^{OZ{(l) + (1-12) 

where C plays the role of a fluctuation and of a background field. 
If we define Sp{(j)) as the bare action and Zp{(p) = e~ ''^'^' , we have: 



Z = I dpCf,{<P)Zpi(t)) = ifio * ■■■ * (Mp-1 * (Mp * Zp) (1.13) 

Let us define 

Zj = {fij* ...*{fip^i*{fip* Zp) ; Zj^i = Hj^i* Zj] (1.14) 

and 

5,-(0) = - log Z, ((/>). 

We see that constructing the ultraviolet limit is the same as finding a bare action Sp{(j)) such 
that the effective, or renormalized action Sq{(J)) converges as /9 ^ oo. 

Remark that the C^ satisfy approximate scaling laws: C^^^{x) ~ M'^C^{Mx). They 
would even satisfy exact scaling if we consider a massless theory with m = which is of- 
ten done. In this case we can use the scaling C^^^{x) = C^{Mx), to perform the famous 
trivial but slightly confusing rescalings of the renormalization group. Defining (J)m{x) = 
M^'^~'^^''^(j){Mx), and Zj^cp) = Zj{(j)j^j) we obtain indeed the equation 

Zj{4>) = (mi * Zj+i){4)M-^) 

This defines the {j independent) TZ operation on the action S as the composition of 4 steps: 
-spatial rescaling by M~^ 



-field rescaling by M-('^-2)/2 

-integration over a fixed (scale 1) fluctuation slice 

-taking the logarithm of the result to reexpress it as an action for the next step. 

In this way the problem of the ultra-violet limit reduces to convergence of the p times 
iterated transformation TZoTZ...oTZSp as p — > oo. The problem of iteration of a fixed map can 
be analyzed as a discrete dynamical system: if a fixed point appears, for instance the free 
field or a theory close to it, one can hopefully trust a perturbative analysis of the vicinity of 
this fixed point. 

In particular we see that for instance polynomial terms such as j (f)'^{x)d'^x scale with a 
factor (corresponding to power counting) which is M^^""" after j steps. For n > 4 they are 
"irrelevant" and die out in the RG flow; the (/>^ term is marginal, and its flow is governed 
by the sign of the "bubble graph" , hence by the sign of the beta function at small coupling. 
Finally the mass term, quadratically growing, is "relevant" . Adding derivative couplings, we 
find that only the wave function term J{V(j))'^{x)d'^x is not irrelevant, but marginal. From this 
very simple analysis follows the classification of renormalizable theories into asymptotically 
free (Gaussian fixed points) and not asymptotically free, and the discovery of non-Gaussian 
fixed points of the RG close to Gaussians, if one modifies for instance slightly the canonical 
scaling of the fields (mimicking non-integer dimensions) for an asymptotically free model. 

Study of infrared critical points by RG is very similar to that of ultraviolet limits, except 
that in an infrared problem, the ultraviolet cutoff is fixed, hence it is natural to give it an 
index or 1; and the slices run towards momentum p = 0, so that as j grows, it indexes 
smaller and smaller momenta. 

Sometimes here a confusion arises: the renormalization group does not simply "exchange" 
infrared and ultraviolet limits. The basic fact to keep in mind is that it flows always in the 
same direction: from small spatial scales to larger ones. Indeed it deduces macroscopic actions 
from microscopic ones, as is traditional in the physical analysis of a phenomenon, where large 
scale effects are explained by smaller ones Q. Hence the direction of flow of RG never changes, 
only the point of view of the observer. In an ultraviolet problem in some sense it is the source 
of the RG river which goes to inflnity, whereas in an infrared problem it is the mouth of the 
river which goes to infinity relatively to the observer. 

1.5 Constructive RG is necessary! 

After this blitz overview of the standard renormalization group of textbooks remark that 
although it clarifies enormously the ultraviolet problem, it is not yet formulated in a correct 
mathematical way. In particular starting form any polynomial action it creates an effective 
action which is obviously no longer polynomial, and this even after a single step! Therefore 
the large field problem (integration on cp, or Q, appears! More precisely the behavior of 
Sgff{(f)) at large (p is unclear, so that starting from a stable interaction, even the second step 
of the RG may be already ill-defined. This point has to be stressed to physicists! 

The solution is to reconsider the single step of the RG, hence the theory in a finite 
slice with a background field and to analyze it in a way which can be iterated correctly 

^We are aware that this traditional view is put in question by duahty in string theory, but in this lecture 
we nevertheless stick to it! 



mathematically. This is the "constructive" version of the RG, as introduced in the 80's 
by Gallavotti and coworkers BCGNOPS, and developed in [GKl-2]. It has been now well 
developed and systematized into a coherent mathematical formalism, by Pordt, Brydges 
and collaborators [P][B]. A related formalism is the "multiscale phase space expansion" of 
[FMRSl-2], now redesigned into a more transparent and explicit formalism in [A][AR2]. 
The main difference is that in the phase space expansion formalism the rescalings are not 
performed, so that the iteration of the RG steps is fully developed. Phase space expansion 
therefore leads to an explicit description of the solution of the RG induction, somewhat like 
Zimmermann's forest formulas "solves" Bogoliubov's induction. 

In these formalisms the main tool consists in performing the key step of integration 
of fluctuation (the mathematically well defined version being called a "cluster expansion") 
and taking the logarithm (the mathematically well defined version being called a "Mayer 
expansion") only in regions where the background field is small. In the other so-called large 
field regions one simply "waits and sees" . Although a full presentation of the constructive 
RG formalism for bosonic theories is beyond the scope of this course, we give in the next 
section in detail the treatment of a single scale model with a cluster + Mayer expansion. 

The constructive RG will be an iteration of this treatment, with the added complication 
(alas very nontrivial) that one cannot expand the "large field regions" and must exploit the 
fact that their weights are small in the probabilistic sense. 

Using the constructive RG approach or the related "multiscale" expansions some models 
have been built, giving the first concrete examples of renormalizable quantum field theories 
fulfilling all Wightman axioms, such as the Gross-Neveu2 model [GKl] [FMRSl]. New models 
not perturbatively renormalizable but asymptotically safe become also accessible, such as the 
Gross-Neveu3 model built by P. da Veiga and collaborators [dV]. In the infrared regime 
bosonic models of renormalizable power counting such as the critical (massless) 01 with an 
infrared cutoff [GK2] [FMRS2], or 4 dimensional weakly self avoiding polymers have been 
controlled [IMl], and their asymptotics at large distance established. Nonperturbative mass 
generation has been established in the Gross-Neveu2 model and in the nonlinear a model at 
large number of components [KMR][K]. Finally the RG when applied to condensed matter 
give rise to many rigorous results and programs, and this is the subject of the last section. 

2 Single Scale Bosonic Model: the cluster and Mayer expan- 
sions 

2.1 Tree and Forest formulas 

In statistical mechanics, the key step of the thermodynamic limit is to take a logarithm to 
pass to intensive quantities which correspond to connected quantities. The minimal discrete 
connecting structures between points are trees. It is therefore worth devoting some time to 
the precise combinatorics of trees, and to the way to generate them systematically between 
points. 

Let n > 2 be an integer, /„ = {l,...,n}, Vn = {{i,j}/i,j £ In,i ^ j} (the set of 
unordered pairs in /„). An element I of Vn will be called a link, a subset of "P^, a graph. A 



graph ^ = {li, . . . ,It-} containing no loops, i.e. no subset {{11,12}, {^2, ^3}, ■ ■ ■ , {ik,h}} with 
k > 2 elements, is called a forest. A connected forest is called a tree. A forest is a union of 
disconnected trees 1, the supports of which are disjoint subsets of /„ called the connected 
components or clusters of ^. 

The first basic result on trees goes back to the XlXth century: 

Theorem 1 (Cayley's Theorem) There are exactly n"~^ labeled trees on a set of n points. 

Exercise 2.1: Prove Cayley's theorem by induction, using the multinomial expansion, and 
showing that there are n!/ Yl{di — 1)! trees with coordination number di at vertex i, using the 
fact that trees always have "leaves". X 

Let S be the space of smooth functions of n{n — l)/2 variables x = {xi)^^^ associated to 
the links of "P^. A Forest formula is a Taylor formula with integral remainder, which expands 
a function H oi S to search for the explicit dependence on given link variables x;. Taylor 
formulas with remainders in general are provided with a "stopping rule" and forest formulas 
stop at the level of connected sets. This means that two points which are joined by a link are 
treated as a single block. (More sophisticated formulas with higher stopping rules are very 
useful for higher particle irreducibility analysis, or renormalization group computations but 
are no longer forest formulas in the strict sense [AR2]). 

Any forest is a union of connected trees. Therefore any forest formula has an associated 
tree formula for its connected components. And therefore, at least formally, any forest formula 
solves the problem of computing normalized correlation functions. Indeed applying the forest 
formula to the functional integral for the unnormalized functions, the connected functions 
are simply given by the connected pieces of the forest formula, hence by the corresponding 
tree formula. It is in this sense that forest formulas exactly solve the well known snag that 
makes connected functions difficult to compute. This snag is that since typically there are 
many trees in a graph, one does not know "which one to choose" when one tries to compute 
connected functions in the (desirable) form of tree sums. Any forest formula gives a particular 
answer to that problem, associating a weakening factor w to the links which remain underived 
(the potential "loops"). This weakening factor w tells us exactly by how much our pondered 
"tree choice" has "weakened" the remaining loop lines. 

Several forest formulas exist, with different rules for w; they correspond to different ways 
of letting forests grow. In one logic, the most "symmetric" one, the forest grows in a random 
way: it leads to a formula first established by Brydges and Kennedy [BK]. In another logic, 
which leads to the "rooted formula" each tree grows layer by layer from a preferred root 
[ARl]. Our presentation here follows [ARl]. 

Applied to an element H oi S, the Taylor Rooted Forest formula takes the form (the 
vector with all entries equal to 1 being denoted by 1): 

Symmetric Taylor Forest Formula 



where the summation extends over all possible forests S, including r = hence the empty 
forest. To each link of ^ is attached a variable of integration wi. The vector ^^(w) is the 
value at which we evaluate the derivative of H; it is the vector {xi)^^^ defined by xi = wf{w), 
where the weakening factor w^--,{w) is mf{wi,l G L^{ij}} if L^{ij} is the unique path in the 
forest ^ connecting i to j, and is if no such path exists (hence if i and j belong to different 
clusters of^. 

The notation u-forest stands for regular "unordered" forests. Decomposing the integration 
domain into r! subdomains according to a complete ordering of the parameters Wi we obtain 
a related formula where the sum runs over so called "ordered" forests in which the links are 
ordered, and the integration parameters Wi follow the ordering. 

Exercise 2.2 Write down the general ordered forest formula. Remark that for n = 2 it 
reduces to the ordinary Taylor formula H{1) = H{0) + /q xH'{wx)dw. Performing a series 
of ordinary Taylor interpolations at each step restricted to the links joining different existing 
clusters, prove by induction the ordered formula, from which the unordered one follows easily 
by regluing the integration domains for the w parameters. Jft 

Exercise 2.3 Write down the ordered and unordered formulas for n = 3 and 4 (very instruc- 
tive). Jl» 

Although we won't use it let us stress that these formulas are not unique. For instance 
there exists an absolutely identical formula, the rooted formula, with a different rule for the 
weakening factor w, now called ^t;f^^(w) (the superscript R standing for "rooted") . It is a 
less symmetric formula since we have to give a rule for choosing a root in each cluster. For 
each non empty subset or cluster C of /„, choose re, the least element in the natural ordering 
of /„ = {1, . . . , n}, to be the root of all the trees with support C that appear in the following 
expansion. Now if i is in some tree 1 with support C we call the height of i the number of 
links in the unique path of the tree 1 that goes from i to the root re- We denote it by /^(i). 
The set of points i with a fixed height k is called the k-th. layer of the tree. The Rooted 
Taylor Forest Formula is then absolutely identical to the Symmetric one, except changing 
the superscript S to R (as "rooted") and defining the weakening parameter w differently, 
by the following rule: 

Rooted weakening factors 

It; ' -I (w) = if i and j are not connected by the S". If i and j fall in the support C of the 
same tree 1 of ^ then 

^t;f^-,(w) = if \l'^{i) — l'^{j)\ > 2 (i and j in distant layers) 

^t;f^-,(w) = 1 if l'^{i) = l'^{j) (i and j in the same layer) 

Wf^-Avf) = W{jj/} if P{i) — 1 = l'^{j) = l'^{i'), ind {ii'} G 1. (i and j in neighboring 
layers, i' is then unique and is called the ancestor of i in %). In particular, if {ij} G S", then 

Exercise 2.4 Prove the rooted formula. Jft 

It remains to see which formulas can be used in the constructive sense and for which 
theories. A very interesting property of the symmetric formula is that it preserves positivity 
in a certain sense: 



Theorem 2 (Positivity of the Symmetric Formula) 

Extended to a symmetric m,atrix by the convention wf^{v^) = IVi, the matrix wfAw) is 
positive. 

Exercise 2.5 Prove this theorem (hint: show that wfA'w) is a convex combination of block 
matrices with 1 everywhere). X 

2.2 Cluster Expansion 

Let us apply the previous formula to study the thermodynamic limit of the single slice cj)^ 
theory. The idea of the cluster expansion is that since perturbation theory diverges we 
must keep most of it in the form of functional integrals. However we can test whether distant 
regions of space are joined or not by propagators, and this will allow to rewrite the theory as a 
polymer gas (with hardcore interactions). When the coupling constant is small, the activities 
for the polymers are small, and the technique of the Mayer expansion which compares the 
hardcore gas to a perfect gas, allows to perform the thermodynamic limit. 

For instance consider the free bosonic (massless) Gaussian measure dfxc ^^ 21 defined by 
the single slice covariance 

Cix,y)= [' ^ ^e-l^-^l^/^- (2.4) 

This propagator is obviously integrable in y at fixed x. 

Put now the regular interaction e~ J^ ^^'dx in a finite volume A and for simplicity let 
us study the pressure 

p= lim^logZ(A) , (2.5) 

where we recall that the partition function Z{K) in a finite volume A is 

Z{K) = /d^c(0)e"^^A'^'(-)'^- . (2.6) 



Let us explain how the Taylor formula (2.3) performs the task of rewriting the partition 
function as a dilute gas of clusters with hard core interaction. 

The set /„ is defined as a partition of A into (unit size) cubes, and clusters are subsets 
of such cubes. We write A = Ui^i^bi, where each 6j is a unit cube, and define Xb as the 
characteristic function of b, and XA = X^ie/ /tbr Since the interaction lies entirely within A, 
the covariance C can be replaced by Ca = XAix)C{x,y)xA{y) without changing the value of 
Z(A). Moreover Ca can be interpolated, defining for I = {i,j} £ Vn 



n 

CAi{xi)i^j:,J{x,y) = J2xb,ix)C{x,y)xb,iy) 

i=l 



+ Yl X{ij}iXb,{x)C{x,y)xb,{y)+Xb,{x)C{x,y)xb,{y)) (2.6) 

{ij}G-Pn 

Remark that Ca{1) = Ca. Now we apply the Taylor formula (2.3) with the function H being 
the partition function obtained by replacing in (2.6) the covariance C by C\{{xi)i^.p ). Here 
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it is crucial to use the positivity theorem, in order for the interpolated covariance to remain 
positive, hence for the corresponding normalized Gaussian measure to remain well defined. 

From the rules of Gaussian integration of polynomials [G J] , we can compute the effect of 
deriving with respect to a given xi parameter, and we obtain that (2.3) in this case takes the 
form 

Z(A) = H{t) = J2 /rf/"CA(x,(w))(0)(n I'dwi) 

where 6j and bj are the two ends of the line I. Since both the local interaction and the 
covariance as a matrix factorize over the clusters of the forest ^, the corresponding contribu- 
tions in (2.7) themselves factorize, which means that (2.7) can also be rewritten as a gas of 
non-overlapping clusters, each of which has an amplitude given by a tree formula: 

Z{A)= fdf^cMe''^^^= ^ flA{Y,) (2.8) 

•^ sets {Yi,...,Yn} j = l 

^(^)= E ill f' dwi) f df,cAX.(^^m 

% on Y l<=%-^'^ ■' 

{ n /''-''«'.W^^'.WC(.,.)j^^}e-/.*'«- (2.9) 

where the sum is over trees X which connect together the set Y , hence have exactly |y| — 1 
elements (if \Y\ = 1, X = connects Y). The measure d^c'y({Xi{h)})('/') is the normalized 
Gaussian measure with (positive) covariance 

CY{X^{^)){x,y) = XY{x){w^{v^){x,y))C{x,y)xY{y) (2.10) 

where w%{'w){x^ y) is 1 if a; and y belong to the same cube, and otherwise it is the infimum of 
the parameters wi for I in the unique path Li{b{x), b{y)) which in the tree 1 joins the cube 
b{x) containing x to the cube b{y) containing y. 

Exercise 2.6 Prove that if Ojj and bij are two positive matrices, their Hadamard product 
Cij = aijbij is again a positive matrix (hint: use square roots) Complete the proof that 
CY{Xi{w)){x,y) is, as announced, a positive covariance, hence that it has a well defined 
associated Gaussian measure. X 

The following bound now summarizes that polymer activities are small enough so that 
one can sum over all polymers containing a fixed point (and absorb a fixed constant per cube 
of the polymer). This will be used in the next section. 

Theorem 3 (Bound on Polymer Activities) 

Given any constant K , for small enough A with Re A > we have 

J2 \A{Y)\k\''\ < 1 (2.11) 

Y such that OeY 
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Exercise 2.7 Prove this theorem (hint: use Formula (2.9). One needs notations to compute 
the action of the functional derivatives in (2.9) by Leibniz rule, the result being cumbersome. 
Then it is useful to remark that the propagators exponential decay (see (2.4)) can absorb 
any "local factorial" of the coordination numbers di of the tree T. This allows to achieve the 
proof 0. Jl» 

2.3 The Mayer expansion 

The Mayer expansion allows to deduce from (2.11) the existence and e.g. the Borel summa- 
bility in A of thermodynamic functions such as the pressure p. 

In the cluster expansion (2.8), the condition that the disjoint union of all clusters is A is 
a global annoying constraint. Remark that the polymer amplitudes are translation invariant. 
In particular the trivial amplitude of a singleton cluster Y = {b} is a number ^o independent 
of b. Redefining Ar{Y) = A{Y)/A^P and Zr{A) = Z(A)/yl|f' we quotient out all the trivial 
clusters so that 

n 

Zr{A) = l + J2 E UMYi) (2.12) 

n>l sets {yj Yn} i=l 

This is the partition function of a polymer gas: the sums over individual polymers would 
be independent were it not for the hard core constraints Yi CiYj = 0. Adding in an infinite 
number of vanishing terms, we can replace the sum in (2.12) by a sum over ordered sequences 
(Yi, ...,Yn) of polymers with hard core interaction and a symmetrizing factor 1/n! coming 
from the replacement of sets by sequences. 

1 " 

Z,(A) = 1 + ^- Y. UMYi) n ^(YuYj) (2.13) 

n>l * sequences (Yi,...,Yn) 2 = 1 l!^^<j5:^ 

\Yi\>2 

where the two-body hard core interaction is expressed by the factors ri{X, y ) = 1 if XnY = 0, 
and ??(X, Y) = a X nY ^ (/>. To factorize these hardcore interactions we apply again the 
symmetric Taylor forest formula (2.3)! More precisely for a fixed sequence (Yi, . . . ,Yn) of 
polymers, we define I„ as the set of these polymers, and define ei = e^^jj = r]{Yi, Yj) — 1, for 
i ^ j- We consider the function 

leVn 
so that H{1) = Ui<i<j<nViyi,yj)- Rewrite (2.13) as 

1 " 

72>1 * sequences (Yi,...,Yn) i=l 

\Yi\>2 



^This method is identical to the one in [R], part III. But remark that although the full amplitudes A{Y) 
defined in (2.9) must be identical to those in [R], the subcontributions associated to particular trees are 
different, since the formula used in [R] was not the symmetric one. 
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1 + E ;^ E ft MY.) E (n f dn^) (n ^0 (n(i + -f (w)^. 

n>l ■ sequences (yi,...,y„) i = l J ^ig;? '^ ^ ^€5 ^ ^i^S 

|y,i>2 
= E;^fEn E (UMYi))c^in,...,Y,)y (2.15) 



n>0 ' fc>l ' sequences (Yi,...,Yi.) i=l 

iyii>2 



where 



c^(ii,...,n) 



1 tree on 



E (uf'dwM n ^0 n (i+^f(wh) (2.16) 



where, as before, wf{'w) is, if / = {ij}, the infimum of the parameters wi' for /' in the unique 
path L^{ij} which in the tree 1 joins i to j. 

Exercise 2.8 Check that the connecting factor C [Yi, ...,Yk) does not depend of the partic- 
ular tree formula, since we have 



C^{Y^,...,Yk)= Yl U'l (2-17) 



G connected graph [f^Cj! 
on{l,....,fc} 



We obtain immediately that 

logZr{A) = J2^ E (i[AriY,))c^iY^,...,Yk) (2.18) 



k 



kl 

fc>l sequences (Yi,...,Yf^) ^^=l 



These formulas can be used together with (2.11) to control the thermodynamic limit 
p = Ao + lim jr^d m log Zr{A) (hint [R]: Every tree coefficient forces the necessary overlaps 
along the links of the tree, and is bounded by 1, since |(1 + ■wf{w)ei)\ < 1. Start from the 
di — 1 leaves of the tree hooked to a given polymer of index i, fix the di — 1 squares of overlap 
of these leaves with the polymer (at a cost lYil'^'^^), then sum over the leaves using the decay 
of (2.11) in the size of polymers. Use the factor i^._i\\ in Cayley's theorem (Exercise 2.1) 

and the fact that X^d \Yi\'^^^^ / {di — 1)! = e'^*'. Then absorb the factor el^*l in the constant 
K of (2.11). Finally iterate to conclude, fixing the last sum to contain the origin, since this 
cancels out the volume factor |A| in (2.5), up to boundary terms which vanish as |A| -^ IR .). 

Exercise 2.9 Prove Borel summability of the theory in a single slice, using the Nevanlinna 
theorem of P.da Veiga's lectures (hint: the analyticity for ReA > and A small is easy. For 
the K^nl bounds on the Taylor remainder at order n, simply rewrite the same analysis than 
above but with n explicit vertices expanded first. Jf> 

A full renormalizable model such as the infrared critical point of ^l can be built by 
extending this single slice analysis to the full multislice model. Let us sketch very briefiy how 
this model would be treated in the multislice expansion of [AR2]. The analysis decomposes 
into three main steps. 
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The multislice model contains propagators C^ in each momentum shce, and space has to 
be therefore decomposed at each shce into a scaled lattice (M^'Z)'^. The union Uj{M^'Z)^ is 
the full phase space for the theory, each cube being some kind of degree of freedom for the 
appropriate field (p^ . The first expansion step is a "small versus large field decomposition" 
which tests, for each cube b G Llj{M^'Z)* of scale j{b), whether a quantity such as Jhiip^^')^ 
is small or not. The maximal connected regions of large field cubes will be treated as single 
blocks, i.e. no cluster expansion will be performed inside them. This step, although perhaps 
not strictly necessary, simplifies a lot the bounds later. 

The second step is the multiscale "cluster expansion" . It has to decouple now all the small 
field cubes and all the large field blocks of all slices j = 0, ...p in (1.10) at the same time. This 
means that it derived both propagators C^ in each momentum slice, this time joining cubes 
of scaled lattices (M^'Z)^. But it has also to derive new links which join the cubes of different 
slices at a same point, this time through vertices. These new "vertex links" can join up to 4 
cubes at the same time, since a vertex can hook to at most 4 propagators in different slices. 
Therefore the symmetric Taylor formula has to be properly extended to accommodate such 
new links [AR2]. Also the formulas in the "momentum slice" direction have to be pushed 
further than simple connectedness, since we want to distinguish convergent activities (those 
with more than 5 external lower momentum legs) from divergent ones (with 2 and 4 external 
lower momentum legs), which require renormalization. As a result the formula in the slice (or 
"vertical" ) direction is not a tree formula, but rather some kind of a "5 particle irreducible" 
formula, with an expansion rule which stops the vertical or vertex interpolations only when 5 
or more links are derived. In spite of these complications, at the end of the multiscale cluster 
expansion the theory is factorized as a gas of polymers swimming in the full phase space, 
hence in Uj{M^'Zi)^, with hardcore interactions. 

The third step is the generalized Mayer expansion, whose role is no longer only to factorize 
and quotient out the normalization (vacuum graphs), but also to factorize and renormalize 
the 2 and 4 point functions. Indeed only the convergent polymers, those which do not con- 
tain subpolymers with 2 and 4 external lower momentum legs, of the previous expansion 
have small summable activities. Divergent polymers require renormalization, the appropri- 
ate counterterms being absorbed into effective constants. This generates more complicated 
formulas, but the logic remains the same (interpolate the hardcore links). The result is some 
complicated but totally explicit formula for expanding the model. It can be viewed as the 
correct constructive analog of Zimmermann's forest formula for perturbative renormalization 
[A]. 

This method, although still very heavy and very technical, is nevertheless clearly more 
explicit and conceptually more transparent than previous methods in which cluster and Mayer 
expansion steps were mixed together in an inductive way (such as in [FMRS2]). 

3 Fermionic theories 

3.1 Introduction 

The initial constructions of renormalizable Fermionic models such as the Gross- Neveu2 model 
[GKl] or [FMRSl] used the same heavy apparatus as Bosonic models, namely slicings in mo- 
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mentum space and direct space, plus cluster and Mayer expansions. As a consequence the 
renormalization group equations obtained were discrete difference equations instead of differ- 
ential equations, a detail which was considered annoying by some physicists. On the other 
hand perturbation theory for Fermion systems is often said to converge, whereas for boson 
systems it is said to diverge. But what does this mean exactly? Unnormalized Fermionic 
perturbation series with cutoffs are not only convergent but entire, whereas Bosonic pertur- 
bation series, even unnormalized and with cutoffs, have zero radius of convergence. Using 
this property there ought to be constructive versions of renormalization for Fermions very 
close to the perturbative concepts. In particular it was advocated recently in [S] that there 
ought to exist continuous flows and differential equations of the renormalization group for 
Fermions. 

Progressively we realized that one can avoid the use of the cluster and Mayer expansions of 
the previous section for Fermionic models (for an early example see [FMRTl]). In particular 
the symmetric forest formula of the previous section, when the interpolated parameters are 
directly applied to the Feynman lines of a graph, gives a particularly simple "three lines" 
construction of interacting Fermions in a single momentum slice [AR3], which is given below. 
In [DR] we applied similar ideas to give a new construction of the Gross-Neveu2 model. This 
construction gives a solution of the theory in the form of explicit sums of finite dimensional 
integrals, containing, however, effective constants which are defined as the (non explicit) 
solutions of the differential renormalization group equations. 

The positivity property of the symmetric Taylor forest formula (Theorem 3 of the previous 
section) was crucial for Bosonic models. But here again it is useful, since it allows to apply 
Gram's estimates: 

Lemma 1 (Interpolated Gram inequality) Let A = aij =< fi ■ gj > be a Gram matrix. 

Gram's inequality remains true for the matrix B = bij = wfJw) < fi ■ gj > , namely: 

IdetBlKllUWllWg.W Vw (3.4) 



Proof: Indeed for fixed ^ and (w) we can take the symmetric square root v of the positive 
matrix w^{w), so that Wij = J2k'^ik''Jkj- Defining the components of the vectors / and g in 
an orthonormal basis for the scalar product < • > to be //" or g^, we define the tensorized 
vectors Fi and Gj with components Fj = Vikfi^ and G^™ = Vjug^ and we have for the 
tensor scalar product < • >t'- bij =< Fi • Gj >t, so that detB < JJi WFiWrHj \\Gj\\t- But 
obviously ||F,||2. = EkmiFt""? = Ekn.vUfr? = wuErnifr? = UW' (we recaU that by 
definition wf^{y/) = 1). □ 

Let us now apply the Symmetric Taylor Forest formula to the computation of the con- 
nected functions of a Fermionic theory. The corresponding Grassmann integral is: 

I J dficii^, V5)P(V^a, V'a)e^('^"-'^") (3.5) 
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where C is the covariance or propagator, P is a particular monomial (set of external fields) 
and S is some general action. We take as simplest example the massive Gross-Neveu model 
with cutoff, for which the action is local and quartic in a certain number of Fermionic fields. 
In a finite box A the action is 

Sa = ^ f dxiJ2Mx)Mx)f (3.6) 

where a runs over some finite set of N "color" indices. A is the coupling constant. There 
may be also a spinor index which we forget in our notations since it plays no role in what 
follows. The covariance C is massive and has an ultraviolet cutoff, hence in Fourier space it 
is for instance r]{p)/{^ + rn) where 77 is a cutoff function on large momenta. We only need 
to know that C is diagonal in color space C(x, a; y, a') = for a 7^ a', and that it can be 
decomposed as C{x,y) = J-^d D{x,t)E{t,y)dt with 

\Cix,y)\<K-—^ (3.7a) 

{l + \x-y\)P 

f \D{x,t)\^dH<K ■ f \E{x,t)\'^dH<K (3.7b) 

for some constant K. 

This decomposition amounts roughly to defining square roots of the covariance in momen- 
tum space and prove their square integrability. It is usually easy for any reasonable cutoff 
model. For instance if r/ is a positive function, we can define D in Fourier space as , J 3^1/4 

ana ii as (^2+^2)3/4 • 

It is the Fermionic covariance or propagator C =< ipa{x)'4>a{y) > which is interpolated 
with the forest formulas. We obtain, for instance for the pressure that the formal power 
series in the coupling constant is a sum over trees on {1, ....,n}, with a distinguished vertex, 
1, sitting at the origin to break translation invariance (similar formulas with external fields 
of course hold for the connected functions): 



3.2 Fermionic Tree Expansion 

p= hm T^LogZiA) = lim -L(^ f d^^c{iJ,i^)e''^^^-'^-^) 
A— >oo |A| A-+C10 |A| J 

^(A7iV«n!)^^e(T,a)(n / dwi) / dxi...dx„<5(xi = 0) 



n=0 1,0- c /ei' 

ra+1 



n C^(^i){xi,yj) i / dfiCx^(^)i'>P^^) n '>Pixi(r),a{r))'>P{xj(^r)a'{r)) (3.8) 

l={ij}e% ^ r=l 

where the sum is over trees 1 which connect together the n points xi, ...,Xn- These trees also 
contain an arrow information a which for each line determines which end was a field and 
which was an antifield. The coloring C contains the color indices of each line of the tree and of 
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each remaining field or antifield. This information completely determines the remaining set 
of uncontracted fields Jlri i^{xi{r),a(r))'^i^ j{r),a' (r)) j that is it determines the functions i{r) 
and j{r) which tell to which vertex the remaining fields are hooked, and the coloring index C 
completely determines the functions a{r), a'{r) (hence their color). Finally e{1,a) is some 
sign ±1 which we do not need to compute. 

Exercise 3.1 Find first a similar but simpler formula for the number of connected graphs 
of the (^^ theory at order n as a sum over trees with weakening factors. (Counting graphs 
with their correct combinatoric factor is equivalent to field theory in zero dimension). Check 
explicitly in the cases n = 3 and n = 4 how the weakening factors when integrated restore 
the right combinatorics for the ordinary graphs (very instructive). Jf> 

Exercise 3.2 Check that (3.8) is indeed the result of applying the Symmetric Taylor Forest 
formula to the Grassmann functional integral (3.5), interpolating the Fermionic covariance 
C and using the Grassmannian rules of "integration by parts". For the courageous reader: 
find an explicit formula for the sign e(T, a). X 

Lemma 2 The index C in the sum (3.8) runs over a set of exactly 2'^N'^^^ elements. 

Proof: At each vertex the circulation of color indices is fixed by a factor 2 (which tells 
which of the two fields -0 is paired with a ip i.e. forced to have the same color by the form 
(3.6) of the action. Once these circulation rules are fixed, the attribution of color indices 
costs N"^ for the first vertex, and by induction a factor N for each vertex of the tree. Indeed 
climbing inductively into the tree layer by layer, at every vertex there is one color already 
fixed by the line joining the vertex to the root, hence one remaining color to fix, except for 
the root, for which two colors have to be fixed. This proves the lemma. D 

Then the expression 

n+l 
C^MCx^(w) (^' ^) n i^iXi{r),a{r))^(^Jir),a'ir)) (3-9) 

is nothing but an n + 1 by n + 1 determinant, with matrix element the covariance 

wf{w)6{a{r),a'{r'))C{xi(^r),Xj(^r')) (3.10) 

between the line r and the column r', where I = {i{r)^j{r')). 

Let us now use these Fermionic tree formulas for proofs of convergence. 

3.3 Convergence of the tree formulas 

A typical constructive result for this Gross-Neveu model with cutoff is to prove: 

Theorem 4 The pressure and the connected functions of the cut-off Gross-Neveu model are 
analytic in \ in a disk of radius R independent of N . 
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Proof: Apply Gram's inequality to the loop lines in (3.8). By Lemma I and (3.7b), this 
determinant is bounded by K'^^^. The spatial integrals are bounded, using (3.7a) by K'^~^. 
The number of trees on n vertices is n"~^ by Cayley's theorem. In the sum over colors, the 
coloring factor A^"+^ of Lemma 2 almost cancels with the factor 1/N^ in (3.8). Hence the 
power series in A for the pressure is bounded by A^ "^^, K'^ for some constant K'. O 

This is perhaps the shortest and most transparent proof of constructive theory yetR The 
theorem above is interesting not only for the analysis of the Gross-Neveu model but also for 
that of the two-dimensional interacting Fermions considered in the next section. In this latter 
case, the "colors" correspond to angular sectors on the Fermi sphere and the factor 1/A^ in 
the coupling is provided by power counting. 

3.4 Renormalization, an overview 

In this subsection we give a brief summary of how to remove the ultraviolet cutoff rj and 
perform renormalization of the Gross-Neveu2 model, following [DR]. 

The covariance C in Fourier space with infrared cutoff A (this is no longer the volume!) 
and ultraviolet cutoff Aq can be written as 



p^ + m^ J A 







Our goal is to prove that the vertex functions of the theory have a non perturbative 
limit as Aq — > cxd, and that they are the Borel sum of their renormalized power series in the 
renormalized coupling constant (because of renormalons they cannot be analytic). 

The main idea is to apply the interpolation forest formula not directly to the lines of the 
Feynman graphs, but to the continuous scale parameter a introduced in (3.11). In this way 
an ordered forest formula is built, in which the ordering of the tree lines plays the role of 
the necessary ordering of momenta in any phase space analysis, namely to distinguish higher 
momenta from lower momenta. This expansion alone however, does not put explicitly into 
display the divergent subgraphs with two and four external legs. An additional construction 
(classes of "chains" in [DR]) has to be performed, to expand further the remaining loop 
determinant in (3.8) so that this structure becomes visible. This remaining expansion has 
to be performed with some care. Indeed we cannot simply sum over the attributions of 
the a parameters of the loop lines to the bands defined by the intervals left between the 
a parameters of the tree lines. This would make visible not only the 2 and 4 point high 
energy subgraphs but also all the other high energy subgraphs, with arbitrary number of 
external legs. And it is a general rule of this kind of mathematical problems that information 
has a price. This particular information would lead to uncontrolled divergences at large 
order. The solution found in [DR] is to reglue together many band attributions into so-called 

^One can also prove the same result using the rooted formula rather than the symmetric one [AR3]. 
Although the proof is somewhat longer, it is interesting since the weakening factor in the rooted formula 
completely factorizes out of the determinant. This second formula may therefore be useful in problems for 
which Gram's inequality is not applicable and the method of "comparison of rows and columns" of [IM2] or 
[FMRTl] has to be used. This method roughly corresponds to Taylor expanding around a middle point further 
and further when fields or antifields accumulate in any given cube of unit size of a lattice covering IR''. 
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"chain classes" so as to expand much less the loop determinant. After this regluing, the only 
information on a high energy subgraph displayed by the expansion is whether this subgraph 
has 2, 4, or more than 4 low energy legs, and this information alone does not introduce large 
order divergences. 

After this subtle point is settled, the rest of the analysis is completely similar to pertur- 
bative renormalization theory. Counterterms for momentum values are introduced which 
subtract the two and four point high energy subgraphs and create flows for the three rele- 
vant or marginal operators of this theory, namely the coupling constant, the wave function 
constant and the mass. Technically the subtractions are performed in direct x space, since 
again we do not want to know the exact structure of loop lines, which would be necessary for 
Zimmermann's subtraction scheme in momentum space. 

Remark that the initial ordering of the tree lines in this formalism reminds of the ordering 
of Hepp's sectors in perturbative renormalization theory (see [R], and references therein); but 
it orders only about "half" the lines of a graph (namely the tree lines) whereas Hepp's sectors 
ordered all the lines of a graph. 

The formalism has also some similarities to the Gallavotti-Nicolo tree expansion [GN], 
but beware that the basic trees in [DR] are not the Gallavotti-Nicolo trees, which correspond 
rather to the auxiliary trees considered for the construction of the "chain classes" . 

4 Renormalization Group in Condensed Matter 

4.1 Many Interacting Fermions and the BCS Problem 

Conducting electrons in a metal at low temperature are well described by Fermi liquid theory. 
However we know that the Fermi liquid theory is not valid until temperature 0. Indeed below 
the BCS critical temperature the dressed electrons or holes which are the excitations of the 
Fermi liquid are bound into Cooper pairs and the metal becomes superconducting. 

During the last ten years a program has been designed to investigate rigorously this 
phenomenon by means of field theoretic methods [BG][FT1-2][FMRT1][S]. In particular the 
renormalization group of Wilson and followers has been extended to models with surface 
singularities such as the Fermi surface. The ultimate goal is to create a mathematically rig- 
orous theory of the BCS transition and of similar phenomena of solid state physics. From 
the start we know that nonperturbative effects must be incorporated in the analysis of su- 
perconductivity since the BCS gap is nonperturbative. But this model is the only one we 
know of in which nonperturbative effects should be accessible in the near future to rigorous 
mathematical control, without any ad hoc modifications (such as introducing an artificially 
large number of colors). Indeed an amazing property of this model is that angular variables 
around the Fermi surface play a role analogous to that of colors, so that an expansion of 
the 1/A^ type should control the BCS regime [FMRT2], in which ordinary perturbation is no 
longer valid. We can call this situation a "dynamical 1/A^" effect. 

Nevertheless the full construction of the BCS model is a long and difficult program which 
requires to glue together several ingredients. The main new idea was to extend the renor- 
malization group of Wison (which analyzes the point singularity p = in momentum space) 
to more general extended singularities. This very natural and general idea is susceptible 
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of many applications in various domains, including field theory (in Minkowski space). But 
this idea has also to be combined with rigorous control of spontaneous continuous symmetry 
breaking, and this generates a lot of technical complications. In these lectures we therefore 
restrict us to the study of the concept of RG around the Fermi surface in a simpler model 
which is at a temperature higher than the temperature of the BCS phase transition. 

The key feature which differentiates electrons in condensed matter from Euclidean field 
theory, and makes the subject in a way mathematically richer, is that Lorentz invariance is 
broken, and density is finite. The field theory formalism remains the best tool to isolate the 
fundamental issues such as the existence of non-perturbative effects. Imaginary (Euclidean) 
time (in the form of a circle, with antiperiodic b.c. for Fermions) corresponds to finite 
temperature. The fundamental state of the theory corresponds to the limit of the temperature 
going to 0. In this limit the imaginary time circle grows to IR Finite density creates the 
Fermi sea and surface. In the simplest case where rotation invariance is preserved (isotropic 
jellium) this surface is simply a sphere. 

The free Fermi liquid theory is therefore defined by Fermion fields with two spin indices, 
and propagator 

Cah{k) = Sah ^ ^-.^ ; e{k) = ^-fi (4.1) 

ikq — e[k) ^"^ 

where a,b G {1,2} are the spin indices. The vector A: is in d spatial dimensions. Adding 
one time dimension there are really d + 1 dimensions. The parameters m and fi correspond 
to the effective mass and the chemical potential (which fixes the Fermi energy). To simplify 
notation we put 2m = 1, /i = 1 so e{k) = fc^ — 1. 

This propagator is rotation invariant, a feature which simplifies considerably the study 
of the renormalization group fiows after branching the interaction. In particular it has a 
spherical Fermi surface. This jellium isotropic model is realistic for instance in solid state 
physics in the limit of weak electron densities (where the Fermi surface becomes approximately 
spherical), but in general a propagator with a more complicated energy function e{k) (such 
as a lattice Laplacian) has to be considered. This is not necessary for our purpose here which 
is to test the constructive validity of perturbation theory. 

Since we work at finite temperature and since Fermionic fields have to satisfy antiperiodic 
boundary conditions, the component /cq can take only discrete values (called the Matsubara 
frequencies) : 

'» = ±^' <^-^) 

so the integral over ko is really a discrete sum over n. As A;o 7^ Vn the denominator in 
C{k) can never be 0. This is why the temperature provides a natural infrared cut-off. When 
T ^ ko becomes a continuous variable and the propagator diverges on the Fermi surface, 
defined by /cq = and \k\ = 1. 

The interaction term in the action is defined by: 

Sj, = ^ f d^xiJ2i^i^f (4.3) 
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Physically this interaction may represent the effective interaction due to phonons. A more 
relaistic interaction would not be completely local to include the short range nature of the 
phonon propagator, but we can consider the local action (4.3) as an idealization which cap- 
tures all essential mathematical difficulties. 

4.2 Renormalization around the Fermi Surface 

It is convenient to add a continuous ultraviolet cut-off (at a fixed scale) to the propagator 
(4.1) because it makes its Fourier transformed kernel in position space well defined, and 
because a non relativistic theory does not make sense anyway at high energies. 

The basic difference between this theory and an ordinary critical point in statistical me- 
chanics or field theory, is that the singularity of the propagator is of codimension 2 in the 
d + 1 dimensional space-time. This changes dramatically the power counting. Instead of 
changing with dimension, power counting in this kind of models is essentially independent 
of the dimension, and is the one of a just renormalizable theory. This can be understood 
basically in the following way. In a graph with 4 external legs, there are n vertices, 2n — 2 
internal lines and L = n — 1 independent loops. Each independent loop momentum gives 
rise to two transverse variables kQ and \k\ and d — 1 inessential bounded angles. Hence the 
2L = 2(n — 1) dimensions of integration for the loop momenta exactly balance the 2n — 2 
singularities of the internal propagators, as is the case in a just renormalizable theory. 

To justify this very crude picture, it turns out that it is very convenient to further de- 
compose the propagator into discrete slices and each slice into discrete angular sectors^: 

C = Y.C^ ■ C^{k) = -l^ (4.4) 

j=i ^ko - e[k) 

where 

f{k) = f[M'^[ki + e{kf)) (4.4) 

effectively forces \ikQ — e{k)\ ~ M~^ . The function / is in Cq^{[1, M^]). The parameter M is 
strictly bigger than one so that the slices pinch more and more the Fermi surface as j — > oo. 
The slice propagator is further decomposed into sectors: 

C(^\k) = Y: C^'^'ik) ; C^'^^ik) = -P^. (4.5) 

<xesi *^0 - e(fc) 

where S-' is a set of roughly M^'^^^'^ angular patches, called sectors, which cover the Fermi 
sphere, with linear dimensions of order M~^ . For instance if d = 2 we simply cut the circle 
into intervals of length 2ttM~^ . 

The RG philosophy applied to this problem is now clear. As before, higher slices give rise 
to local effects relatively to lower slices. After more and more slice integrations one obtains 
effective actions which govern longer and longer distance physics. These effective actions are 
however more complicated than in the field theory context. In rotation invariant models such 

'''A continuous slicing in the style of [DR] is also possible; but the discrete slicing is more in Wilson's spirit. 
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as the one above, renormalization of the two point function can be absorbed in a change 
of normalization of the Fermi radius. It removes all infinities from perturbation theory at 
generic momenta [FTl]. But the flow for the four point function is a flow for an infinite 
set of coupling constants describing the momentum zero channel of the Cooper pairs [FT2] . 
In the case of attractive interaction for A, if the temperature cutoff is lowered to zero, this 
flow diverges at the BCS scale, where the symmetry linked to particle number conservation 
is spontaneously broken, giving rise to the effective BCS theory for the Cooper pairs. 

Like in the previous section the key problem from the constructive point of view is to 
understand the resummation of perturbation theory in a single slice, in order to iterate the 
renormalization group step. 

Curiously, although power counting does not depend on the dimension, momentum con- 
servation in terms of sectors in a fixed slice depends on it and this has crucial constructive 
consequences. In d = 2 we have the "rhombus rule": four sectors meeting at a vertex must 
be roughly two by two equal^. This means that two dimensional condensed matter in a slice 
is directly analogous to an A^-vector model in which angles on the Fermi surface play the 
role of colors [FMRT3]. This allows to complete the proof that the radius of convergence 
of perturbation theory is independent of the slice index j. Roughly speaking, the model in 
slice j is a vector model with A^ = M^ colors. Power counting at a vertex costs M"^^ for the 
space-time integration, and the propagators in fixed sectors earn a scaling factor M^"^^ each. 
Since there are in average two propagators per vertex, we remain with a factor M^^ = 1/N 
left per vertex, which is exactly what is necessary to pay for the color sums and obtain a 
uniform radius of convergence by Theorem 4 of section 3. This completes the sketch of the 
proof of: 

Theorem 5 [FMRTl] In two dimensions, there exists a finite k > independent of 
j such that the power series in A of the Schwinger function for the interacting Fermionic 
measure with propagator (4-5) and interaction (4-3) has a convergence radius of at least k. 

Furthermore it ought to be possible in two dimensions to resum perturbation theory to 
build an effective action until the BCS temperature. As an example, the convergent part 
of the expansion (which does not include 2 or 4 point subgraphs) has been resummed in 
[FMRTl]. 

In three dimensions, two momenta at a vertex in a slice do not determine the third and 
fourth: there is an additional torsion angle, since four momenta of same length adding to 
are not necessarily coplanar. The radius of convergence of perturbation theory is still 
independent of the slice index j, so that Theorem 5 also holds in three dimension [MR], but 
it is much harder to prove than in dimension 2, and until now it is not clear that this partial 
result alone allows a full constructive analysis of the model up to the scale where the BCS 
symmetry breaking takes place. 

In dimension 1 the Fermi surface reduces to two points, and there is also no proper BCS 
theory since there is no continuous symmetry breaking in two dimensions (by the "Mermin- 
Wagner theorem"). Nevertheless the many Fermion system in 1 dimension gives rise to an 

®A precise version is given in [FMRTl]. Remark that there is a logarithmic correction to this rhombus 
rule, due to the case ol a nearly collapsed rhombus. A nice improvement to avoid further problems with this 
logarithm is to define angular sectors longer in the tangential than in the radial direction. 
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interesting non-trivial behavior, that of the Luttinger Uquid [BG]. 

4.3 The Weakly Coupled Anderson Model 

This model describes a single electron in a random potential. This is not strictly speaking a 
field theory. But renormalization group and field theoretic methods can also be applied to 
models such as self repelling walks [IMl]. This model lies in this category. The two-point 
function is given by 

S{x,y;E;XV;e) = — ^— . {x,y) 



1 



( 2 ^ , ,J —^ X -77 (^-y) (4-5) 



where r/ is an ultraviolet cutoff and F is a random potential (multiplication operator in 
x-space), for instance distributed with a Gaussian (regularized) white noise for which the 
covariance is a (regularized) delta function. Indeed in this model the singularity in Fourier 
space lies on the surface \p\ = vE, just as for the interacting Fermions of the previous 
subsection. The perturbative expansion for the averaged Green's function is a resolvent 
expansion in which the integration over the potential creates structures similar to the 4)\ 
graphs, as for self avoiding polymers; the main difference is in the combinatorics for the 
graphs, which is the one of a 'W = component" theory (in particular vacuum parts are 
forbidden). 

The main mathematical question of physical interest is to prove the existence or non- 
existence of localized states. This in turn amounts to study the behavior of the average of 
the modulus square of the two-point function 

Zim,_o j \S{x,y-E-\V-e)\'' d^{V) (4.6) 

The rigorous results in more than one dimension are restricted up to now to the strong 
coupling regime, where it has been proved that this average of the modulus square of the 
two-point function decays exponentially [FS],[AM]. This in turn implies localization of all 
states. 

At small coupling it is expected, that 

-i) in two dimensions the average of the modulus square of the two-point function decays 
exponentially with a rate dependent on the coupling 

- ii) in three dimensions it decays only like a power which means that some states are 
delocalized. 

One expects also that for all (non zero) couplings the mean value of the two-point function 



Urn, 



e-+0 



J S{x,y;E;XV;e)dfi{V) (4.7) 



decreases exponentially, and that it is real analytic in E. 
As a first step in this direction, we have proved [MPR]: 
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Theorem 6 In two dimensions, there exists some k > such that for A small enough 

S{x,y;E;XV; X^^^) dfi{V) ~ e-"^*(") ^' l^'^l (4.8) 



which is the expected rate of decrease . 

To suppress the regularization e = A^"'"'* in our theorem, hence to prove the exponential 
decay of (4.7) is well under way. 

We can only give a flavor of the arguments used in Theorem 6. A key new ingredient is to 
use the rhombus rule in two dimensions for the part of the theory in which the momenta are 
close to the singularity p^ = E, to create loops with approximate momentum conservation. 
These loops are then smaller than expected because of a Ward identity which is somewhat 
involved technically [MPR]. 

But a more immediately accessible and instructive point is to apply the decomposition 
of the critical surface into discrete angular sectors also to this problem. One finds that the 
random potential sandwiched between sector cutoffs of a given slice 

Vaip) V{p - q) r]^'{q) ~ Va,a' for p £ a, q£ a' (4.9) 

is a random matrix whose elements are labeled by the shell cells. In two dimensions, since 
all momenta in a white noise distribution are equiprobable, and since the momentum p — q 
of V approximately determines p and q on the circle \p\ = \fE (because of the rhombus 
rule), this random matrix is approximately a self adjoint matrix with random independently 
distributed entries 

Via, a') = V(cj',ct) = V(-a, -a') • (4.10) 

Therefore we can use the mathematical theory of such matrices for instance to bound the 
deviation or probability that eigenvalues become anomalously large [M] . In three dimensions 
this iV by A^ matrix is more complicated because 

V^^^t = Vry iff a + a' = t + t' (4.11) 

so that approximately yiV matrix elements corresponding to the same momentum transfer 
are equal. The study of such matrices with nonindependent entries, and more generally the 
study of the weakly coupled 3 dimensional Anderson model remains a mathematical challenge. 

5 Further topics 

In this last section we give a brief list of open problems for constructive methods and rigorous 
renormalization group studies. 

- Mass generation in the Gross-Neveu model without cutoff. Prove that at large 
number A^ of colors, the model without ultraviolet cutofjhas a spontaneously generated mass, 
for instance by "gluing together" the analysis of the ultraviolet limit in [DR] and the mass 
generation with ultraviolet cutoff in [KMR]. This could be called "constructive dimensional 
transmutation", after Coleman. 
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- Ultraviolet limit of the nonlinear a model in tw^o dimension. This is a long 
standing problem. Although there are partial results, for instance for the hierarchical model, 
there is no clear rigorous construction using asymptotic freedom of this model. Once it is 
built, there is also the problem of gluing the ultraviolet analysis to mass generation (at large 
N) to obtain "dimensional transmutation" also in this model. This would be interesting in 
view of the controversial issues raised by Patrascioiu and Seller on this kind of model. 

- Fermi liquid and BCS theory in dimension 2: Since [FMRTl], the full construction 
of the BCS vacuum at zero temperature can be considered a sound mathematical constructive 
program, although its technical realization is very hard, and implies a tricky constructive 
analysis of the infrared problems associated with the Golstone boson of the BCS continuous 
symmetry breaking. 

A more accessible task is to precise the mathematical status of Fermi liquid theory it- 
self. This is important also in view of the debate around the nature of high-temperature 
superconductivity. Fermi liquid theory is not valid at zero temperature because of the BCS 
instability. Even when the dominant electron interaction is repulsive, the Kohn-Luttinger 
instabilities prevent the Fermi liquid theory to be generically valid until zero temperature. 
There are nevertheless two proposals for a mathematically rigorous Fermi liquid theory: 

- one can block the BCS and Kohn-Luttinger instabilities by considering models in which 
the Fermi surface is not invariant under p -^ —p (we suggest to call this, for obvious reasons, 
the "egg model"). In two dimensions it should be possible to prove (nonperturbatively) that 
in this case the Fermi liquid theory remains valid at zero temperature, and the corresponding 
program is well under way [FKLT]. This program requires to control rigorously the stability 
of a non-spherical Fermi surface under the renormalization group flow, a difficult technical 
issue [FST]. 

- a simpler proposal, advocated in [S], is to study the Fermi liquid theory at finite tem- 
perature above the BCS transition temperature. As seen above, the temperature acts as an 
infrared cutoff on the propagator in the field theory description of the model. Hence in this 
point of view, the nontrivial theorem consists in showing that stability (i.e. summability of 
perturbation theory) holds for all temperatures higher than a certain critical temperature 
whose dependence in terms of the initial interaction is proved to be of the correct BCS form 
[S]. 

- Fermi liquid and BCS theory in dimension 3: Here the constructive analog of 
[FMRTl] must be found first. [MR] is only a first step in this direction. Having worked for 
several years with J. Magnen on this problem, it seems to me the most beautiful I met in 
constructive theory until now. It is well posed and simple. For instance one could ask: is 
the sum of all convergent contributions to the theory analytic in a finite disk? Perturbative 
power counting suggests that the answer should be "yes" , but because of the torsion angle, 
the constructive problem is surprisingly hard. Again an understanding of this constructive 
issue should clear up the way for many subprograms, such as the Fermi liquid theory in 3 
dimensions for "egg models" or above the BCS transition, and the BCS theory itself. 

Anderson Model, Anderson-Mott phase transition in dimension 3 See the pre- 
vious subsection. 

- Continuation to Minkowski space, constructive scattering theory Here again 
we feel that the key for instance to a constructive analysis of asymptotic completeness in 
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field theory might be related to renormalization group analysis around a surface, this time 
the mass shell. The "constructive return" to Minkowski space is of course full of interest 
and of important issues for the physics of the models (scattering, bound states, time depen- 
dent problems). Constructive study of time dependent problems in statistical mechanics or 
condensed matter theory is also clearly a potential field for many developments. 

- In the future one could focus more systematically on clarifying the mathematical status 
of the fast developed theories of the recent decades (Conformal theories, Non Abelian Gauge 
Theories, in particular Supersymmetric ones. Topological Theories). Some constructive tools, 
of course presumably with unexpected additions, might be of help. In this line one can also 
dream of Constructive Duality and Constructive String or Membrane Theory. 
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